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Abstract We study magnetothermal instability in the 
ionized plasmas including the effects of Ohmic, ambipo- 
lar and Hall diffusion. Magnetic field in the single fluid 
approximation does not allow transverse thermal con- 
densations, however, non-ideal effects highly diminish 
the stabilizing role of the magnetic field in thermally 
unstable plasmas. Therefore, enhanced growth rate of 
thermal condensation modes in the presence of the dif- 
fusion mechanisms speed up the rate of structure for- 
mation. 
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1 Introduction 

Role of the magnetic field in the dynamics of the 
gaseous astrophysical systems is generally studied 
within the framework of ideal MHD equations. In this 
simplified approach, which is good under certain condi- 
tions and circumstances, it is assumed that the coupling 
between the charged and neutral species of the system 
is perfect. But the perfect coupling assumption can be 
violated, in particular, when the density of the charged 
particles can be much lower than that of the neutral 
species. For example, one shoul d note to this fact 



dense molecular clouds (e.g., Ciolek and Robergel 



2002). Some authors have also criticized application of 
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ideal MHD equation in modeling accretion discs, par- 
ticular for th e discs around young stellar objects (e.g., 
WardldEil . 

The very existence of various species with differ- 
ent masses and electrical charges and their collisions 
and the possible momentum transfer between the par- 
ticles should be considered in any theory of struc- 
ture formation in interstellar medium (ISM). How- 
ever, interest towards modeling astrophysical plasmas 
within a mu ltifluid appr o ach has been raised over recent 



years (e.g.. iFalld l2003t lO'Sullivan and Downed 



2007 



Pandev and Wardld 



Pinto Galli. and BacciottiTkoO; 

20081 : Inoue and Inutsukal 120081 ). In particular, simu- 
lating multifluid systems is a challe nging area. For ex- 



ample, IChoi. Kim, and Wiital |2009h incorporating am- 
bipolar diffusion in the strong coupling approximation 
into a multidimensional magnetohydrodynamic (MHD) 
code based on the total variation diminishing scheme. 
More generalized Multifluid numerical schemes arc also 
studied during recent years by some authors (e.g. JFalld 
20031; lO 'Sullivan and Downesll2007l ). 



On the theoretical side, efforts to understand physics 



2007 



of multifluid plasmas are in progres s (e.g-. lTassis and Mous choviasl 
Kunz and Mouschoviasteooih . Recently. Pandev and Wardle| 



(hereafter PW) clarified the relationship be- 
tween the fully ionized and weakly ionized limits by 
developing a unified single-fluid framework for the dy- 
namics of the plasmas of arbitrary ionization. In 
another related study, general expressions for the 
resistivities, the diffusion timcscales and the heat- 
ing rates in a three-fluid medi um are obtained by 
Pinto. Galli. and Bacciotti l|2008t ). They showed that 



the value of the Ohmic resistivity is increased in a col- 
lapsing cloud and the ambipolar diffusion occurs on a 
time scale comparable to the dynamical time scale. 

Among various physical mechanisms responsible for 
the density inhomogencitics in the ISM, it has been re- 
alized for a long time that thermal instability can be 



e fficient processes (e.g., Fieldl 1965 : Burkert and Lin 



2000; 
2008 : 



Hennebelle and Audit 



Gazol. Luis, and Kim 



200 



200 



Inoue and Inutsuka 
Thermal instabil- 
ity in a cooling medium can also be affected when dy- 
namics of di fferent charged species (e.g., du st particles) 
are included ( Ibahez and Shchekinovl2002h . It is shown 
that negatively charged particles stimulate the thermal 
instability in the sense that the conditions for the in- 
stability to hold are wide r than similar conditions in 
a single-fluid description (IKopp and Shchekinovl 12007 : 



Shukla and SandberdlioO 

Thermal instability within ambipolar regim e also 
studied by Neiad-Asghar and Ghanbaril ( 2003 ) and 
Neiad-Asgharl|2007l ). where in the former the frictional 
heating by the ion-neutral is included . In an other sim- 
ilar study, IStiele. Lesch. and Heitschl (|2006l ) revisited 
the problem of clump formation due to thermal instabil- 
ities in a weakly ionized plasma. However, dynamics of 
ions and their contribution to the net cooling function 
are not included in these studies o f thermal instability 



in am bipolar regime. Therefore, iFukue and Kamava 
( 20071 ) extended the linear classical thermal instability 
to a case, in which dynamics of neutrals and ions and 
their interactions are considered. They showed that 
ion-neutral friction with the magnetic field affects the 
morphology and evolution of the interstellar matter. 
Thermal instability for a system ob eying generalized 
Ohm' s law has also been studied by iBora and Talwar 
(|l993f ). They found that the instability criterion in- 
volves the field strength, resistivity and electron inertia 
terms for transverse perturbations, but for the parallel 
to the ambient magnetic field the instability criterion 
is independent of all these non-ideal effects. However, 
they did not apply their analysis to the structure for- 
mation due to the thermal instability in the interstellar 
medium. 

In this paper, our goal is to include non-ideal effects 
in a magnetized thermally unstable plasma. We follow 
a recent approach by PW, in which ions, electrons and 
neutrals are included. Then, non-ideal Ohmic, Hall and 
ambipolar terms appear in a generalized form of the in- 
duction equation. We study role of these terms in mag- 
netothermal instability in the linear regime. The ba- 
sic assumptions and the equations are presented in the 
next section. Linearized equations and the dispersion 
relation are derived in section [3] Analyzing effects of 
non-ideal mechanisms on the thermally unstable modes 
are done in section[4] We conclude by a summary of the 
results and possible implications in the final section. 



2 Basic multifluid equations 

Basic assumptions and the equations of our model are 
constructed based on the approach of PW. Here, an 



energy equation is introduced as well. The possible ef- 
fects of charged dust particles on th e thermal instability 



are n eglected for simplicity (e.g., IShadmehri and Dib 
2009). We assume that the system is consisted of ions, 
electrons and neutrals. However, it is difficult to con- 
sider the full set of dynamical equations for a multifluid 
system with three different species. But one can reduce 
the set of multifluid equations to a more manageable 
set of equations with some extra terms due to the non- 
ideal effects as has been done by some authors (e.g., 
PW). The basic equations of our analysis are discussed 
here, but one can refer to PW for a detailed discussion 
about the assumptions and the key approximations. 

For each component of the system, the continuity 
equation is written as 



dp] 
dt 



V • (pjvj) = 0, 



(1) 



where pj = rrijUj is the mass density, v 3 is the velocity, 
and rij and uij are the number density and particle 
mass of ions, electrons and neutrals, i.e. j = i, e, n. 
The continuity equation for the bulk fluid is obtained by 
summing up equation {l} for each species. Therefore, 



dp 

m 



V ■ (pv) = 0, 



(2) 



where p = Pj ~ Pi + Pn and v = (p,Vi + p n v n )/p are 
the bulk fluid density and velocity, respectively. 

We shall assume that the ions are singly charged and 
adopt charge neutrality, i.e. n\ = n c . The momentum 
equations for the electrons, ions and neutrals are 



-VP e - n e e(E + — x B) 
c 

-p c ^v cj (w c -w ] ) 1 (3) 



(9v - V' 
ft(-5T + Vi ' Vvi ) = ~ VPi + ?li<3 ( E + - x B ) 

at c 



'Pi ^'fa - V J')) ( 4 ) 



J=e,n 

Pn(-^+V,rVv n ) = -VPn+jOn ^ ^(Vj-Vn), (5) 

j=c,i 

respectively. The electron and ion momentum equa- 
tions (J2])-© contain on the right hand side pressure gra- 
dient, Lorentz force and collision momentum exchange 
terms where Pj is the pressure, E and B are the elec- 
tric and magnetic field, c is the speed of light, and i/jk 
is the collision frequency for jth component with fcth 
component (i.e., p 3 v 3k = pkVkj)- 
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The bulk momentum equation can be derived by 
adding equations ^3j)— (JSJ) to obtain 



,dv . p[p n 
p(— + V Vv) + V- ( V D V D ) 

at p 



-VP+-JxB, (6) 

c 



where P = P + Pj + P n is the total pressure, vd = 
Vj— v n is the ion-neutral drift velocity, and J = n e e(vi — 
v e ) is the current density. PW showed that the vdvd 
term in equation ^ can be neglected for dynamical 
frequencies satisfying 



to < 



D(3 C 



y/plpZ^l + Dft 



(7) 



where D = p n / p is the neutral density fraction and (3 C 
is the ratio of the cyclotron frequency of the electron 
oj c to the sum of the electron-ion and electron-neutral 
collision frequency v c (PW), i.e. (3 C = <jJ c /v c or 



(eB/m e c) 



(8) 



Then, we can recover the single-fluid momentum equa- 
tion as 



3v 



Vv) = —VP + - J x B. 



(9) 



To obtain an equation for the evolution of the mag- 
netic field, we need to drive an expression for the elec- 
tric field E in terms of the fluid properties to insert into 
Faraday's law. We present the main steps for obtain- 
ing the final form of the induction equation following 
detailed and extensive calculations of PW. First, we ob- 
tain a relation for vd by rewriting the ion and neutral 
equations of motion ((3J and (j4|) as 



Pencil J V D 



-V(P C + Pi) 



-M^ + vt.VvO 



and 



,<9v n 



(10) 



(pi^in + Pe*4n)vD = - P^(~gf + V n-Vv n ) 



-VP n 



(11) 



Multiplying equation (flO|) by p n and equation (fTT|) by 
Pi and then adding and noting that p c v cn <C PiV m , we 
obtain 

VB = D iJ^ + YIjL. D lL + i ^±. (12) 

cp\v m piv m pifin Pe en e 



Obviously, the above equation is not valid when i> m 
tends to zero unless we relax the assumption that we 
used, i.e. p c v n <C Pi^m- In the weakly ionized limit, 
equation (|12[) reduces to the strong coupling approxi- 
mation, i.e. vd ~ (J x B) / (cpiV- m ) . 

Now, we need a relation for the electric field. The 
electron momentum equation in the zero electron 
inertia limit, yields an expression for the electric field 
in the rest frame of the ions 



E=--xB 

c 



(:?7 C 



cen. 



-vd, (13) 



where a = e 2 n e /m e (i' en +i' e i) is the Ohmic conductivity 
and J is given by Ampere's law, J = (c/47r)V x B. But 
we note Vj = v + Pvd- Substituting equation Q12[) 
into equation (|13|) . the induction equation dTi/dt = 
-cV x E becomes (PW) 



?? = Vx[(vx B)-^J-^JxB 

dt LV ' c c 



(JxB)x B], 



(14) 



where B = B/P> and the Ohmic, Hall and ambipolar 
coefficients arc ?y = c 2 /47rcr, t^h = cP/47ren G and r/A = 
D 2 B 2 /AirpiVin, respectively. 

Now, if we write the energy equation for ions, elec- 
trons and neutrals and add up the equations, and ne- 
glecting the external cooling and heating of electrons, 
we obtain the bulk energy equation that is given by its 
usual form as follows 



7 - 1 v dt 



VP) 



7 



-PV 



- P n-V.[KuV\\T + K_ L W±T}=0, 



(15) 



where O is the net cooling function. Also, the coefficient 
of thermal conductivity K has the values K || and K± 
in directions parallel and perpendicular to the magnetic 
field B. 

There are some points regarding the energy equa- 
tion (|15p and the net cooling function, in particular. 
In analysing thermal instabili ty, the net cooling func- 
tion, Q, has a vital role (e.g., Fieldl 1965 ). In addition 
to the usual heating-cooling terms, three extra heat- 
ing terms appear in the net cooling function because 
of the non-ideal effects, notably Ohmic, Hall and Am- 
bipolar. While the usual heating-cooling terms are ac- 
tually local functions and depend only on the physical 
quantities of the system, non-ideal heating terms are 
determined mainly by the spatial variations of the mag- 
netic field. In other words, all non-ideal heating terms 
vanish for a system with uniform magnetic field. In 
such circumstances, equilibrium states of the system arc 



•1 



solely determined by the local heating-cooling terms. It 
means that the classical criteria of thermal instability 
is still applicable when the magnetic field is uniform. 
But growth rate of the thermally unstable modes may 
significantly modify due to the non-ideal effects on the 
distribution of magnetic field through our modified in- 
duction equation. On the other hand, linear analysis of 
thermal instability is restricted to keeping terms which 
are of the first order and since non-ideal heating terms 
are of the second order, they do not appear in the lin- 
earized set of the equations, irrespective of the profile 
of the magnetic field. Thus, possible effects of non- 
ideal effects on the growth rate of thermally unstable 
modes appears through modified induction equation in 
our simplified analysis which is restricted to the linear 
regime with uniform initial conditions. 

Therefore, equations ©, ©, (O and ((TSJ) along 
with the ideal gas equation of state, p = (R/ p)pT, are 
the main equations of our model in order to study ther- 
mal instability in a multifluid system. The equations 
look like similar to the ideal MHD equations, except for 
the induction equation where some extra terms emerge 
due to the non-ideal effects such as ion-electron or ion- 
neutral drift velocities. 



Pa Pa T 



(21) 



where the derivative Cl p = (dil/dp)T and = 
(dil/dT) p are evaluated for the equilibrium state. 

We introduce the coordinate system e x , e y , and e 2 
specified by 



B B x k 

e,-— , e y - - , e x -e y xe z . {22) 



Also, we introduce the following wavenumbers 
k P = Kl - l)po^p(i?c s T ) _1 , 
fc T = /i(7-l)^T(i?c s )- 1 , 

= - l)-K"||] -1 (-RcsPo), 



(23) 



Now, we can write the dispersion equation using the 
following non-dimensional quantities, 



Tp = T' (XT = T> <TK " = k^ l 



k 



(24) 



3 Linear perturbations 

In the local homogeneous equilibrium state, we have 
p = pa, P = Pa, T = T , B = B , v = 0, and 
ft(po,To) = 0. We assume perturbations of the form 



X( r > t) = xi exp(wi + ik ■ r), 



(16) 



where xi is the amplitude of the perturbations, uj is the 
growth rate of the perturbations and k is the wavenum- 
ber of the perturbations. We are interested in modes 
where uj is real and positive, i.e. condensation modes. 
Then the linearized equations arc 



Our stud y differs fro m the classical thermal instability 
analysis (|Fieldlll965l) in introducing a generalized form 
for the induction equation because of the non-ideal ef- 
fects. For numerical purposes, it is more convenient, 
to re-write Ohmic, ambipolar and Hall coefficients in 
terms of some non-dimensional parameters. In doing 
so, we have 



° s A 



(25) 



(26) 



ujpi + ip k ■ vi = 0, 



(17) 



(27) 



uj Po v 1 + zkPi + ^Bo ■ Bi)— - i(k • B )5i = 0, (18) 

47T 47T 

-Pi - t 7t — Pi + Pa^pPi + PattrTi 

7-1 (7 - 1 )po 



where the non-dimensional parameters 0, A and Ti are 
defined as 



O = 



kpC s 



P e (m e /m*)uj c 



(28) 



+{K\\k\ + K 1 .k\)T 1 = 0, 



(19) 



wBi + iB (k • vi) - i(k ■ B )x 1 - r?k x (k x Bi) 
+r ?A kx{[(kxB 1 )xB ]xB } 

-7 ?H (k-B )(kxB 1 ) = 0, (20) 



(29) 
(30) 
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where a = (va/c s ) 2 and i>a is the Alfven velocity. Here, 
the effective ion mass is m* = p/n c . Then, the Hall fre- 
quency is defined as uh = eB/m*c. Also, the cyclotron 
frequency of electrons is written as lo cc = eB/m e c. 
Therefore, the characteristic equation becomes 

Y 7 + P 6 Y 6 + P 5 Y 5 + PiY A + P 3 Y 3 + P 2 Y 2 

+P 1 Y + P = 0, (31) 

where Y = uj/k p c s and the coefficients are 

a 2 e 



Po 



7 



■(ctt + cr K - <7 P 



7 



2a^O + ^-(o T + o K - o p )H 2 + a£(a 



2. 

1 



X(o T + K ) Op 



P 3 = U 2 + [(1 + OO + (2a£ + ^^)(o T + o K ) 

7 

-(— VpM + O 2 + [K + a + -){o T + o K ) 



7 



7 



P 5 =^A 2 + (1 + 0(0 + o T + o K )A + O 2 
+2(o T + o K )C + iU 2 + a + at + 1 



7 



(32) 



Pi = —(o T + ok - o p )[20 + (£ + 1)4] + (33) 

7 



£ 1 + £ 

P 2 = -(0- T + CTK - O p )A 2 + [( -)(o T + K 

7 7 

-o p )0 + <(£ + 1)]4 +-(o T + ok - o p )0 2 

7 



(34) 



7 

-^]0 + £W 2 +«£(« + 2) (35) 

7 



P 4 = e(^T + ok)A 2 + [(1 + i)(o T + o K )0 + 2ai 
+C + l]A + (or + o K )0 2 + (2 + a + a£)0 

+i(o T + o K )H 2 + (a + a£+ -)(o T + o K ) - — (36) 



(37) 



P 6 = (l+0A + 2O + o T + o K . (38) 

and £ = cos 2 8 and 9 is the angle between Bo and k. 
Also, we have ok = (J K lt ^ + ok ± (1 — C)- O ur analysis 
is based on dispersion equation (f3Tj) and its roots arc 
analyzed in the next section. 



4 Analysis 

If we neglect the non-ideal effects and set O = A = 
H = 0, the algebraic equation (f3"Tj) reduces to the stan- 
dard magne othermal instability characteristic equation 
dFieldlll965l) . It is hard to do Hurwitz analysis in or- 
der to study all possible roots of the dispersion equa- 
tion qualitatively due to its complicated coefficients. 
But equation (|3Tj) is of odd degree in Y and must 
therefore admit at least one positive real root for non- 
perpendicular perturbation if the last term Po is nega- 
tive. Obviously, a positive root implies monotonic in- 
stability, i.e. condensation mode. Therefore, as long as 
different species doe not contribute to the net cooling 
function the condition of instability is the same as the 
ideal case, i.e. Po < or (ot + ok — o p ) < 0. 

Also, we found that equation (p?T|) can not be factor- 
ized to a product of lower order polynomials unless the 
effect of Hall diffusion is neglected. If we set Tt = 0, the 
equation becomes a fifth order polynomial multiply by 
a second order polynomial, i.e. Y 2 + (O + A£,)Y + a£. 
Then, the fifth order polynomial describes unstable 
modes and two stable waves modes are described by 
setting this second order polynomial equal to zero. 

Now, we find real and positive root of equation ((3"T|) 
numerically. We take the parameters 7 = 5/3, a = 1, 
ot/o p = 1/2 and o p ok = 0.01 for comparison to 
Field (1965). Figure Q] shows growth rate versus the 
wavenumber of the perturbations for Ohmic, ambipo- 
lar and Hall regimes, separately. For easier comparison, 
growth rates in the non-magnetized and ideal magne- 
tized cases are shown by dotted and dashed curves, 
respectively. In the top plots of Figure [T] we set 
A = H = and the non-dimensional Ohmic coeffi- 
cient O varies from 0.1 to 5. Each curve is labeled by 
its corresponding O. Growth rate of the condensation 
mode increases with the coefficient O and the effect is 
stronger for the transverse perturbations. Obviously, 
larger growth rate implies a more magnetothermally 
unstable system. So, Ohmic diffusion has a destabiliz- 
ing effect on the condensation mode. As Ohmic coeffi- 
cient becomes large, the profile of the growth rate tends 
to the non-magnetized case. In other words, the effect 
of magnetic field on the condensation mode diminishes 
in the presence of an efficient Ohmic diffusion. 

Middle plots of Figure [T] show growth rate in am- 
bipolar regime, i.e. O = H = but A ^ 0. Here, non- 
dimensional ambipolar diffusion coefficient A varies 
from 0.1 to 10. As in Ohmic regime, a destabilizing ef- 
fect is seen due to ambipolar dissipation. Bottom plots 
of Figure [T] show growth rate of the condensation mode 
when the system is in the Hall regime, i.e. O = A = 
and TC ^ 0. Destabilizing effect of Hall diffusion on 
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Fig. 1 Growth rates of magneto-thermal condensation mode versus wavenumber of the perturbations for three cases: 
Ohmic regime (top), Ambipolar regime (middle) and Hall regime (bottom). Our fiducial parameter are a = 1, 7 = 5/3, 
ot/o"p = 1/2 and a p ax = 0.01. Dotted and dashed curves are representing growth rates in the non-magnetized and ideal 
magnetized cases. Each curve in the top, middle and bottom plots is labeled by fiducial values of O, A and TC, respectively. 
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the thermal instability is amplified with the coefficient 
Ti. The wavelength of perturbation that maximizes the 
growth rate profile shifts to longer wavelengths as non- 
ideal effects become stronger. Then, larger clouds are 
the most likely structures to be formed due to mag- 
netothcrmal instability in comparison to the ideal case 
when non-ideal effects are taken into account. 

Our results are explained using a few simple points. 
Different species not only may contribute to the net 
cooling of the system but their dynamical roles will ap- 
pear either by direct interactions and momentum ex- 
changes or through magnetic field forces for the charged 
particles. However, the bulk density and velocity arc 
mainly determined by the most massive particle, i.e. 
neutrals. Because electrons are very light and their dy- 
namical roles can always be neglected. Moreover, in the 
weakly ionized case, the mass fraction of ions are small 
in comparison to the total density of the thermal gas. 
In the ideal MHD description, all the charged particles 
are well coupled to the magnetic field lines and the neu- 
trals are moving along with the charged species. But 
in reality, magnetic coupling of the charged particles 
may vary depending on the mass of the particles and 
the magnetic strength. Obviously, electrons are more 
coupled to the magnetic field lines in comparison to 
the other charged particles. Thus, the strength of the 
magnetic force may change. In our model, response of 
the system to the perturbations is mainly determined 
by the distribution of the magnetic field lines. So, non- 
ideal effects lead to a reduction to the magnetic pressure 
and this makes the system more unstable. 

Considering the above point, growth rates in Figure 
[1] are physically understandable. Growth rates corre- 
sponding to the non-magnetized and ideal magnetized 
cases are also shown in all panels of this Figure. We 
can consider two thermal and magnetized systems: one 
including non-ideal magnetic effects and the other one 
in the ideal magnetized case. Let's consider transverse 
perturbations, first. If we compare ideally magnetized 
case with non-magnetized situation, we see that mag- 
netic field stabilizes the system, i.e. growth rates de- 
creases in the ideally magnetized case. It can be simply 
explained by the magnetic pressure that provides ad- 
ditional support against perturbations. Thus, in the 
ideally magnetized case, there are actually two types of 
pressures: thermal and magnetic. So, in order to have 
thermal instability, perturbations should overcome to- 
tal pressure which consists of the thermal and the mag- 
netic pressures. Now, we consider non-ideal effects such 
as resistivity or ambipolar. Obviously, these mecha- 
nisms are dissipativc. It means that magnetic flux and 
eventually magnetic pressure is reduced due to the non- 
ideal effects. Thus, there is less total pressure in non- 
ideal case in comparison to the ideal case. This implies 



that growth rates are increased with non-ideal parame- 
ters and the curves are shifting upwards, i.e. from ideal 
case to the non-magnetized case. 

Equation (f2"9")l shows that the ambipolar coefficient is 
directly proportional to the ratio D = p n / p. So, small 
values of A correspond to when D tends to zero. On the 
other hand, one can simply show that O is directly pro- 
portional to D according to equation . In a highly 
ionized case, i.e. D tends to zero, unlike Ohmic and 
ambipolar mechanisms, Hall diffusion does not disap- 
pear in the induction equation. So, the only effect of 
D is seen directly via parameters O and A. Figure [T] 
shows an extensive study of the parameter space, and 
we can see simply that as A or O tends to zero the 
growth rates tend to the classical result of the thermal 
instability. 

We can now apply our linear multifluid magne- 
tothermal instability to structure formation in ISM. As 
an example, we consider structure formation in typ- 

In our calcula- 



ical HI clouds (jWolfire et al 
tion, we adopt m\ = 30to f 



19951). 

and m n 



2.33m p for 

the ion and mean neutral mass, respectively. Here, 
m p = 1.67 x 1CP 24 g is the proton mass. We assume 
that the temperature of HI cloud is 100 K. The number 
density of neutral component and the magnetic field 
are estimated to be 71.9 cm~ 3 and 10~ 6 G. Also, we 
have 7 = 5/3. Under these conditions , CII cooling 



is a d ominant mechanism according to IWolfire et 



1 199,51 ) and the cooling rate depends on the ionization 
degree and the fraction of CI in CII (i.e., /en)- For 
example, for ionization degree 10 -6 and feu = 0.01, 



the cooling rate becomes k p c s = 9.28 x 10 



-14 



(e-i 



Fukue and Kamavall2OO70 . Assuming that all ions are 
singly ionized, our diffusive non-dimensional parame- 
ters becomes O/a = 3.4 x 10~ 10 , A/a = 9.4 and 
H/a — 0.35. Thus, Ohmic dissipation is not operating 
in HI regions. But ambipolar and Hall diffusions are 
significant processes to be considered for HI structure 
formation due to magnetothcrmal instability. 



5 Summary and Conclusions 

We studied magnetothermal instability with the ef- 
fect of non-ideal Ohmic, ambipolar and Hall diffusion. 
Our linear analysis shows that the criteria of insta- 
bility does not change comparing to the ideal case as 
long as charged species do not contribute to the net 
cooling function. Also, the system becomes more un- 
stable in the presence of non-ideal effects and it is 
more probable to have larger clouds comparing to the 
ideal case. Although the vital role of magnetic field 
dissipation in very dense interstellar clouds is a key 
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process in standard theories of star formation (e.g. 



Nakano and Umebavashil 1198 6). our results show that 



such non-ideal mechanisms may operate in thermally 
u nstable systems such a s HI r egions or warm ISM. 
Inoue and Inutsukal (|2008l ) studied formation of 



structures in a weakly ionized and magnetized interstel- 
lar medium using two- fluid magnetohydrodynamic sim- 
ulations. When orientation of magnetic field is perpen- 
dicular to the flow, the rate of form ation of clouds slows 
down significantly according to llnoue and Inutsuka 
(|2008l ). Actually, linear ideal magnetothermal insta- 
bility shows that magnetic field pre vents struct ure for- 
mation transverse to the field lines ( Field 19651 ). How- 
ever, our results show that the stabilizing effect of mag- 
netic field drastically diminishes because of the dissi- 
pative processes like Ohmic or ambipolar. Therefore, 
we think, the opportunity of fast molecular cloud for- 
mation directly from the warm neutral medium would 
highly increase due to the non-ideal effects, at least in 
the linear regime. 

One should note that enhanced growth rate due to 
the non-ideal terms is independent of the true mech- 
anisms of the operating processes. But we think this 
independency is valid in the linear regime and the non- 
linear evolution of the system will depend on the type 
of dominant diffusion process. In particular, Hall dif- 
fusion significantly differs from Ohmic and ambipolar 
diffusions regarding to the energy considerations. In 
fact, it is known that Ohmic and ambipolar diffusion 
are dissipativc processes and reduce magnetic energy of 
the system. But Hall diffusion does not contribute to 
the dissipation of the magnetic energy and so, it is not 
a dissipativc mechanism. Its main role goes back to re- 
distributing the current within the system and it may 
lead to enhanced dissipation of the magnetic energy 
because of Ohmic or ambipolar diffusions. Numerical 
simulations show that the dissipation rate of the MHD 
turbulence is stro ngly affected by the strengt h of am- 
bipolar diffusion ( Choi. Kim, and Wiita 20091 ). But it 
is an open question to explore possible role of Hall diffu- 
sion in no nlinear evolut ion of thermally unstable system 
(see also. lWardlell2004l) . 
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